Statistical properties of single-file diffusion front 
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Statistical properties of the front of a semi-infinite system of single-file diffusion (one dimensional system 
where particles cannot pass each other, but in-between collisions each one independently follow diffusive mo- 
tion) are investigated. Exact as well as asymptotic results are provided for the probability density function of (a) 
the front-position, (b) the maximum of the front-positions, and (c) the first-passage time to a given position. The 
asymptotic laws for the front-position and the maximum front-position are found to be governed by the Fisher- 
Tippett-Gumbel extreme value statistics. The asymptotic properties of the first-passage time is dominated by a 
stretched-exponential tail in the distribution. The farness of the front with the rest of the system is investigated 
by considering (i) the gap from the front to the closest particle, and (ii) the density profile with respect to the 
front-position, and analytical results are provided for late time behaviors. 



Diffusion of a set of labeled particles in the one- 
dimensional (ID) line with the restriction that the particles 
cannot pass one another (so that the order of particle la- 
bels remains same over time) is referred as single-file diffu- 
sion (SFD). The notion of SFD was originally introduced by 
Hodgkin and Keynes to describe the transport of ions through 
narrow channels in biological membranes [1]. Subsequently it 
has been realized in a variety of systems including superionic 
or organic conductors [2] and certain zeolites with a ID pore 
structure 021]. Interacting ID particle systems, which include 
the SFD, are of general interest 0, BE, and have found ap- 
plications in diverse fields ranging from modeling collective 
dynamics of molecular motors [6] to highway traffic flow |7|]. 

The restriction on the motion of particles in SFD makes 
the movement of individual particles (tagged-particles) in the 
bulk anomalously slow. In an infinite system, the probability 
density function (PDF) po(x,t) of finding a particle at posi- 
tion x at time t , initially starting at the origin, approaches the 
Gaussian propagator at long time, 



Po(x,t) ■ 
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exp 
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where the mean-square displacement (x 2 (t)) is proportional 
to \fi, rather than t as found in the case of a simple diffu- 
sion process. This result was first obtained by Harris |Ht], and 
subsequently reestablished by various others using different 
methods fl 111 El EI EI El El El ■ 

The anomalous mean-square displacement in SFD sys- 
tems has been observed experimentally in diffusion measure- 
ments made with pulse field gradient NMR on methane or 
CF4 confined in uni-dimensional channels in AIPO4-5 ze- 
olite fni EM EH- Recent advances in experimental tech- 
niques have made it possible to create well-defined SFD 
model-system by confining paramagnetic colloidal spheres in 
trenches fabricated by photolithography, where the particle- 
particle interaction can be precisely adjusted by an external 
field ifeoll . or by confining charged colloidal particles in ID 
channels generated by scanning optical tweezers [21]. These 
systems allow the trajectories of individual particles to be fol- 
lowed over long period of time with video microscopy, and 



hence it has been possible to experimentally verify the theo- 
retical result given by Eq. (Q}, together with the anomalous 
mean-square displacement 12011 . 

The theoretical and experimental studies in SFD systems, 
thus far, have concerned with the properties of tagged- 
particles in the bulk. However, a particle at the boundary 
should obey quite different laws from those of a tagged- 
particle in the bulk. The theoretical effort to understand those 
laws has been missing thus far, which may be associated with 
the lack of an corresponding ideal experimentally accessible 
model-system until recently. However, with the present level 
of experimental sophistication |20l 21], it would be feasible 
to have accurate measurements on a particle at the boundary 
of a SFD model-system. This motivates the present Letter 
to initiate the needed theoretical study on the properties of 
the front of a semi-infinite SFD system. One indeed finds 
that the SFD-front obeys rather different and interesting laws, 
which would encourage experimental measurements of new 
kind. In a classical single-file system, namely, the Jepsen 
gas ll22tl — the particles undergo binary elastic collisions, but 
move ballistically in-between collisions — the front-velocity 
has been studied recently [23], which is also related to the 
fitness in a simple toy model of biological evolution of quasis- 
pecies 124112511 . Therefore, the theoretical results of this Letter 
on the SFD-front would also be of interest to this discipline. 

Consider a semi-infinite system consists of hard-point 
Brownian particles, initially at t = distributing on the semi- 
infinite real line (—00,0] according to a Poisson process hav- 
ing density (intensity) p, with the "zero"-th particle locating 
at the origin, — i.e., the particle-positions are: xo(0) = and 
jCj+i (0) < xi(Q), i = 0, where the inter-particle gaps, 
/, = x, (0) — Xj + \ (0), i > 0, are independent and identically dis- 
tributed (i.i.d.) according to the PDF = p exp(— p/,). As 
this system evolves in time, each particle diffuses indepen- 
dently having the same diffusion constant D, except the hard- 
core repulsion forbids one particle to pass another one when 
two particles meet, i.e., each particle acts as a hard wall for 
the other. Thus, the order of particle labels is preserved over 
time, Xj+i(t) < x;(f),V? > 0,z"> 0. The aim of this Letter is 
to provide a comprehensive study of the front, described by its 
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position Xo(t), as the above SFD system evolves in time. It ob- 
tains exact as well as asymptotic analytical results for the fol- 
lowing experimentally measurable quantities associated with 
the front: (A) the propagator p(x,t), — the PDF of finding the 
front xo(f) at position x at time t, (B) p max (z,t) — the PDF 
of maximum front-position z = max({xo(T)},0 < X <t), (C) 
Pf(t,z) — the PDF of first-passage time T of the front to a 
given position z (z > 0), (D) i//b(e,f) — the PDF of the gap 
e(t) = Xo(t) —Xi(t), between the front and the next particle, 
and (E) the density profile p(£,t) as a function of the distance 
£, £ > 0, measured from the position of the front. 

We begin by explaining the two main steps involve in the 
calculations. The first step relates the SFD problem to a 
non-interacting one where Brownian trajectories Xj(t), j — 
0, . . . can pass each other |8||. In the non-interacting prob- 
lem, whenever two particles pass each other their order in the 
ID line is interchanged, and hence, the order of the particle 
labels is not preserved over time. The ordered SFD trajecto- 
ries {x,(f)} can be constructed by tracing along the parts of 
non-interacting trajectories {Xj(t)}, keeping the order in the 
ID line fixed, — one single trajectory x,(f) in the SFD may 
consists parts of several non-interacting trajectories {Xj(t)}. 
Therefore, under this mapping, the position of the SFD front 
is given by the position of the rightmost particle in the non- 
interacting system, i.e., 



x (t)=wa(X (t),Xi(t),...,X a> (t)), Vr. 



(2) 



In the second step, one considers time evolution ofN+1 non- 
interacting Brownian particles, starting from an initial config- 
uration at t = where one particle is located at the origin and 
N particles are independently and uniformly distributed in the 
interval [— L,Q] in the real line. In this non-interacting sys- 
tem, the required quantities are calculated for the rightmost 
particle, — which is the front Xo(t) in the SFD system, — av- 
eraging over all possible initial positions of N particles in the 
interval [— L,0]. One finally takes the limit N — > °°, L — > °o, 
while holding the particle density p — N/L fixed, to obtain 
the required quantities for the original semi-infinite SFD sys- 
tem. Note that, the probability = (1 -l/L) N , that an 
inter-particle gap is larger than Z, in the limit N — * °°, L — ► °°, 
keepingp —N/Lfixtd, converges to = exp(— p/), which 
is signature of a Poisson process. Having described the main 
steps, we now proceed with the calculations. 

Propagator. — Consider P(x,t) — f^ oa p(xo,t)dxo, the 
probability that front-position Xo(t) is less than x at time t. In 
the non-interacting particle picture, according to Eq. (0, this 
is exactly the probability that at f, all the particle-positions 
{Xj(t)} are less than x. For a single Brownian particle with 
the initial position X, (0) = yu the probability that at t the posi- 
tion Xj(t) is less thanx, is Pi (x,t\yi) = 1 — <p(Xi~yi,t)dXi, 
where 



(j)(x,t) = 



1 



VAnDt 



exp 



xr 
4D~t 



(3) 



is the single-particle Gaussian propagator. Explicitly, 



P 1 (x,t\y i ) = l 



(4) 



Now averaging over the initial position of y; in [-L, 0] gives 
Pn(x,t) = {i/L) f L P{(x,t\yi)dyj. Therefore, the probability 
of such no n- interacting particles at t having positions less 
than x, after averaging over their initial positions in [— L,0], 
is P N , L (x,t) = [PixCV )] ■ Finally the limit N — > <», L — > oo, 
holding the particle density p — N/L fixed, leads to 




N/L=p 



in which, 



S(u) = exp(— m 2 ) — \/7tuerfc(u). 



(5) 



(6) 



The desired probability P(x,t) is, therefore, obtained by mul- 
tiplying Eq. (0 by Pi (x,t |0), for the particle that started from 
the origin. Explicitly, 



P(x,t)=F 

where the scaling function is, 
erfc(«) 



F(u,t) = 



1- 



exp 



pVDt 



S(u] 



(7) 



(8) 



Equation (0, together with Eqs. ([8]) and (0, gives the exact 
properties of the SFD front for all p\Jl)t. The limit p — re- 
covers the result for a single diffusive particle. The propagator 
of the front is given by, p(x,t) = dP(x,t)/dx. For p VDf 3> 1, 
using the asymptotic form S(u) ~ exp(— u 2 )/2u 2 in Eq. (0, 
one finds the limiting form 



F(u,t) 



( u~a(pVDi/2^n) \ 

\ b(pVDi/2^i) 



in which the scaling parameters are, 

a(c) — Vine — lnlnc/2Vlnc 
b{c) = 1/2^0+ ■■■ , 



(9) 



(10a) 
(10b) 



and the scaling function is the Fisher- Tippett-Gumbel (FTG) 
distribution function 126], 



G(x) — exp(— exp(— x)j , —oo<x<°o. 



(11) 



Equivalent limiting result has been obtained by Arratia 111 111 . 
However, we have arrived at the result here more transpar- 
ent way from the exact result Eq. (0. The propagator of 
the front p(xo,t) for py/Dt ^> 1, after the appropriate scaling 
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x = [x /\f^-a{p\/^/2^n)]/b{p\/Di/2^/n), is given 
by the FTG PDF g{x) = exp(-x - exp(-x)) . 

Maximum front-position. — Let Q{z,t) be the probabil- 
ity that the front of the SFD system, with the initial condi- 
tion x (0) = and -oo < x i+l (0) < x,(0) < 0, i = 1,...,°°, 
does not cross the position z {z > 0), up to time t. Clearly, 
g( z ,f) = Prob[max({x (T)},0 < x < t) < z], and the PDF 
of the maximum of the front-positions in the time inter- 
val [0,f] is given by p max (z,t) = dQ(z,t)/dz. In the non- 
interacting picture, Q(z,t) is exactly the probability that none 
of the particle crosses z up to t [cf. Eq. d2j]. Now, the 
probability that a Brownian particle starting from y (y <. 



0) at t = does not cross z (z > 0) up to t is Q27|, I28H . 
Qi{z,t\y) — erf([z — y]/V4Df), known as the survival prob- 
ability in the literature. Averaging over the initial position 
in [-L,0] gives Q x {z,t) = (l/L) J° L Q UL (z,t\y)dy. Then for 
N non-interacting particles with initial positions averaged in 
[— L,Q], the probability that none of them crosses z up to t 

is simply Qn,l{^J) — [6l(z>0] ■ Finally the limit N — ► °°, 
L — ► oo, holding p =N/L fixed, leads to [cf. Eq. ©] 



N.L- 



N/L=p 



exp 




(12) 



Similar expression with z = 0, i.e., 5(0) = 1, has been encoun- 
tered elsewhere, e.g., see Ref. [29]. However, our interest lies 
in the behavior of Q(z,t) for z > 0, which is exactly the prod- 
uct of the above expression and the probability Q\ (z,f |0), that 
the particle starting from the origin does not cross z in the in- 
terval [0,t]: 



Q(z,t) = a 



where the exact scaling function is, 



-.J 



Q.(u,t) — erf(u) exp 



2pVDt 



S(u) 



(13) 



(14) 



For pyDt ^s> 1, Q,(u,t) has the same limiting FTG distribution 
as in Eq. ©, however, now at different values of the scaling 
parameters, namely, 



Q(i*,f) 



u - a(p \fWt I \J%) 
b{p^Dt/^H) 



(15) 



Therefore, for pVDt 1, the PDF of the maximum position 
f>max(z,0, a f ter tne appropriate scaling, is given by the FTG 
PDF g(x) = exp(— x — exp(— x)). 

First-passage time. — The time x at which the SFD front 
crosses a specified position z for the first time, is the first- 
passage time (FPT) of the front to z. The FPT problems 
have wide implications and therefore, are of general inter- 
est 0, 0]. The probability /" p v [x,z)dx, that the FPT to 
Z (z > 0) is greater than f, is exactly the probability that the 
SFD front has not cross z in the interval [0,t], which is given 



by Q(z,t) in Eq. ([T3J. Therefore, the PDF of the FPT to a 
given position z is pf(t,z) = —dQ(z, x) jdx, with proper nor- 
malization Jq pp(x,z)dx = Q(z,0) — 1. It is convenient to 
express the PDF of the FPT in the form 



Pf(t,z) =A(t,z) exp(-T(T,z)), 
in which the two functions A (x,z) and T(t,z) are, 



(16) 



^^^"(vfe)' <17> 



z 2 2 P Vm c f z 



(18) 



The limit p = reproduces the known result for the FPT of 
a single Brownian particle, as expected. For p ^ 0, while 
X — > limit is still governed by the essential singular behav- 
ior of the single particle, at long time the asymptotic behav- 
ior of pf(x,z) is dominated by the stretched-exponential tail 
exp(— 2p\fDx I \f%), in contrast to the t~ 3 / 2 power-law be- 
havior of the single particle. 

Finally, consider the question: How far is the front at a later 
time, away from the rest of the system? There is, of course, 
no unique way to measure it quantitatively. The remaining 
part of this Letter addresses the question by considering two 
measurable quantities, namely, the gap between the front and 
the next particle e(f) = Xo(t) — X\{t), and the density profile 
with respect to the front. 

Gap. — Consider ¥o(e,t) = /" Vo(£' ,t)de' , the probabil- 
ity that at time t the gap is greater than £. Let R(e,t\x)dx 
be the conditional probability that the gap is greater than £, 
given that the front xq (t ) lies within [x + £ , x + £ + dx] . Clearly, 
Wo(e,t) = /"oo R(e,t \ x)dx. Let us use the notation R^.l (e , f |jc) 
to refer R(s,t\x), in the finite non-interacting system, where 
initially one particle is located at the origin, and others uni- 
formly distributed in [— L, 0] . Then, in the non-interacting pic- 
ture, Rn ,t\x) dx is given by the probability that at t, one 
particle is located in [x + E,x + e + dx], and rest particles 
have positions less than x. The particle in [x + e,x + e + dx] at 
t is either the one initially located at the origin or one out of 
particles initially uniformly distributed in [— L, 0] — the coiTe- 
sponding PDF is q NX (x + £ , t ) = (N/L) J ° L (x + £ - y; , t) dy t . 
Therefore, R NtL (e,t\x) = <j)(x + e,t)P N , L (x,t) + qN,L{x + 
e,t)Pi(x,t\0)Pff-i^(x,t), In the limit N — * oo, L — > oo, holding 
p =N/L fixed, qN,L(x,t) — > (p /2) erfc (x / y/4Dt) . Thus, tak- 
ing this thermodynamic limit inRN.L(s,t\x) yields, R(e,t\x) = 
W(e/V4Di,x/V4Dt,t)/V4Dt, where [cf. Eqs. ©-© and 
©], 



(19) 



W(S,u,t)=-^= exp[ -(u + 8f-^-^S{u) 



pVDt erfc(« + S)F(u,t). 



The desired probability for the gap can be obtained 
by integrating R(e,t\x) over x, — i.e., *Fo(£,i) — 
J-coW(e/V4Dt,u,t) du. For p^/Dt > 1, the main con- 
tribution comes from the second term of Eq. ( fT9l ). Using 
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the limiting form of F(u,t) given by Eq. © and the 
asymptotic form of pVDf erfc(x), one finally obtains 
*Fo(£,f) ~ exp(— p (f)e), for p\J~Di > 1. Here p is the 
effective density near the front, given by 



Po( f ) 



ln(pV5f/20r) 



lV2 



(20) 



Density profile. — The density of particles at a position 
x(t), is the PDF of finding a particle at position x at time t, 
is (p /2) erf c (x A/ 4Df ) + 0(x,f). However, the goal here is 
to compute the density ~p(£,t) at a distance measured from 
the front-position Xo(t), — i.e., the PDF of finding a particle 
at a distance I away from the front-position. Let ~p(£,t\x) 
be the conditional PDF of finding a particle at position x — £, 
£>0, given that the front is at position x. Clearly, ~p(£,t) — 
J-ce~p(£,t\x)dx. In the non-interacting picture, ~p(£,t\x) is 
the joint PDF of finding one particle at position x — £, an- 
other one at x and rest of them having positions less than x. 
Taking the thermodynamic limit in the non-interacting pic- 
ture, one gets, ~p(l,t\x) = qu m (x — £,t) qn m (x,t) P(x,t) H , 

where q\i m {x,t) = (p /2) erfc (x/\/4Dt) , is the PDF of find- 
ing a particle at position x at time t starting from the interval 
(—oo,0), and P(x,t) is given by Eq. (0. The remaining terms 
in p (£, t \x) are the contributions of the events where the trajec- 
tories of the particle starting from the origin reach x and x — £ 
respectively, and can be neglected for large p\fWt. Therefore, 
integrating p(i,t\x) overx, for pVDt ^> 1, 



p(£,t) 



erfc 



h(u,t) du 1 



(21) 



where h(u,t) = p\/Dte.rfc(u)F(u,t). In the large pyDt limit, 
h(u,t), in terms of the scaled variable x = \u — a(cj] / b(c) 
with c = py/Dt jl\f%, converges to the FTG PDF, g(x) = 
exp(— x — exp(— x)). Therefore, in the limit p\fT)t — > as 
b(p*jDt /2^fn\ — > 0, h(u,t) becomes Dirac-delta function at 
u = a(p \f~Dt j r L\f%) , and hence, Eq. (f2Tb yields 



,t) w — erfc 

' ; 2 



U \ 2y/n J y/4D~t 



(22) 



Analogous limiting form also arises for the density of near- 
extreme events in i.i.d. random variables J30ll . In the £ — > oo 
limit, ~p{°°,t) = p, as expected. For £ — 0, using the asymp- 
totic form erfc(a) ~ exp(— a 2 )/ay/n, one recovers Eq. (f20b 

withp(0,f)=P (0- 

In conclusion, this Letter has provided a comprehensive 
study of statistical properties of the front of a semi-infinite 



single-file diffusive system. It would be interesting to investi- 
gate the effect of bias, finite -ranged inter-particle interactions, 
and the disorder presence in the system. The ultimate interest, 
however, lies in realizing the results in experimental systems. 
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